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Abstract 

The structure of l/2 ± states in 13 C up to around the 3a + n threshold (E x = 12.3 MeV) is 
investigated with a full four-body 3a+ n orthogonality condition model (OCM) calculation, where 
the 3a OCM, the model space of which is the subspace of the 3 a + n model, describes well the 
structure of the low-lying states of 12 C including the 2f, Oj, and O4" states, which have been 
recently observed above the Hoyle state (0</)- A full spectrum up to the 1/2/ (l/2j/) state is 
reproduced consistently with the lowest five 1/2“ (three l/2 + ) states of experimental spectrum. It 
is shown that the 1/2./ and 1/2/ states are characterized by the dominant cluster configurations 
of 9 Be(3/2“,l/2“)-|-a, while the ground state 1/2/ has a shell-model-like structure. The observed 
monopole transition strengths to the 1/2/ 3 states are consistently reproduced for the first time. 
These results indicate that the excited 1/2“ states have cluster structures. They are compared to 
those by the previous work with the shell model. On the other hand, the 1/2+ state is found to 
have a loosely bound neutron structure in which the extra neutron moves around 12 C(g.s) core 
with 15 orbit, reflecting the fact that this state appears by 1.9 MeV just below the 12 C(g.s)+n 
threshold, while the 1/2^ and l/2j/ states are characterized by 9 Be+a structures. We found 
that the l/2g" state located above the 3a + n threshold is the Hoyle analogue state in 13 C, the 
wave function of which is described by product states of constituent clusters, (OS') 3 (S) n , with the 
probability of 52 %. 

PACS numbers: 21.10.Dr, 21.10.Gv, 21.60.Gx, 03.75.Hh 
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I. INTRODUCTION 


Cluster as well as mean field pictures are important to understand the structure of light 
nuclei [1-41. A lot of cluster states are known to exist in light nuclei ( 3 ] as well as neutron 


rich nuclei J 5 , 6] and hypernuclei [7j. The typical cluster state is the Hoyle state, the second 
0 + state (O 2 ") at E x = 7.65 MeV in 12 C. This state is located just above the 3a disintegrated 
threshold, and is characterized by the large monopole transition rate sharing about 16 % of 
the en ergy-w eighted sum rule 8]. The microscopic and semi-microscopic cluster models in 
1970’s |9Mlll] demonstrated that the Hoyle state has a loosely coupled 3a cluster structure. In 
2000’s, however, it was found that the Hoyle state has a remarkable aspect of the a-particle 


0 - 



eventually confirmed the second 2 + state ( 2 ^) of 12 C, which had been predicted at a few 
MeV above the Hoyle state by the microscopic cluster models in 1970’s. A new observation 
of the 4 + state at 13.3 MeV was reported in Ref. 26]- In addition, two broad 0 + states, O 3 
and O 4 , were observed 22] at 9.04 MeV and 10.56 MeV with the widths of 1.45 MeV and 
1.42 MeV, respectively, the results of which are consistent with the theoretical predictions 
by the semi-microscopic cluster model, 3a OCM (orthogonality condition model), using 
the complex scaling method jEl, 28]. The underlying structure of these newly observed 
states has been also discussed by the rigid 3a cluster picture 29S, [30|] and the container 
picture 31]. According to the experimental and theoretical analyses 22 , 0 El], the 0 ^ 
state has a prominent 8 Be(0 + )+a structure with a higher nodal behavior, while the Oj 
state is characterized by a linear-chain-like structure having the dominant configuration of 
8 Be(2 + )+a with a relative D -wave motion. 

The discovery of the a-particlc condensate aspect in the Hoyle state has triggered search¬ 
ing for Hoyle-analog states in other light nuclei such as 16 0 32, 33], n B 34-36], 13 C 37, 38], 


and 56 Ni 48] etc. The definition of the Hoyle-analog states is that the constituent clusters 
make nonlocalized motion occupying mainly the lowest orbit (or some excited orbits) of the 
cluster mean-held potential in nucleus. Reflecting their gas-like properties, they seem to 
appear around the disintegrated threshold of the constituent clusters. The structure study 
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of 16 O has recently made a significant pro gres s 


32 


33 


cluster model in 1970’s 


39], following the study of the 12 C+a 


.40 ] and 1980’s 41]. The six lowest 0 + states of 16 0 (including 


the ground state with a doubly closed shell structure) up to E x ~ 16 MeV around the 4 a 


threshold were for the first time reproduced well toge 
decay widths within the framework of the 4a OCM 


her with the monopole strengths and 


32]. The OCM is a semi-microscopic 


cluster model, which is an approximation of RGM (resonating group method) and is exten¬ 


sively described in Ref. 42], Many successful applications of OCM to ordinary nuclei as 
well as hypernuclei are reported in Refs. { 3 , 4, 7]. According to the 4 a OCM calculation, 
the five lowest excited 0 + states were found to have cluster structures with the substantial 
monopole strengths comparable to the single particle monopole strength j^l, 33]. The Ojf 


state at 15.1 MeV just above the 4 a threshold was assigned as a strong candidate of the 4a 
condensate [ill, (OS') 4 , with the probability of 61 %. A novel structure of 12 C(Hoyle)+a 
has been discussed in the highly excited energy region around E x ~ 18 MeV by the po¬ 


tential model 


the 4 a OCM 


44], and the 12 C(0^)+a model 45]. On the other hand, 


the 1/22 state of n B just above the 2 a + t threshold was pointed out to have a Hoyle ana¬ 
logue structure, (OS) 2 (0S)t, with the probability of 60 % 36]. In addition, multi-a gaslike 


states 


12 


46 


47] have been experimentally explored in 36 Ar, 40 Ca and 56 Ni 48]. An 


in¬ 


teresting viewpoint of clustering, a container picture, has quite recently been proposed to 
understand cluster dynamics not only for the gaslike cluster states but also for the clus¬ 
ter states which had been regarded as having localized cluster structures like the inversion 


doublet band in 20 Ne 


49 


50] and linear-chain states 51. 


Isoscalar (IS) monopo 
in the low-energy region 


e transition strengths are very useful to search for cluster states 


33 


34]- The IS monopole excitations to cluster states in light 


nuclei are in general strong as to be comparable with the single particle strength, and their 
experimental strengths share about 20% of the sum rule value in the case of 4 He, n B, 12 C, and 
16 0 etc. js. 33]. As a typical case in light nuclei, the IS monopole strength function of 16 0 52 1 
was discussed up to E x ~ 40 MeV [33j. It was found that 1) two different types of monopole 
excitations exist in 16 O; one is the monopole excitation to cluster states which is dominant in 
the lower energy part, and the other is the monopole excitation of the mean-held type such 
as one-particle one-hole ( lplh ) which is attributed mainly to the higher energy part, and 2) 
this character of the monopole excitations originates from the fact that the ground state of 
16 O with the dominant doubly closed shell structure has a duality of the mean-held-type as 
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well as alpha-clustering character. According to the Bayman-Bohr theorem 53|, the dual 
nature of the ground state seems to be a common feature in light nuclei. The 4a OCM 
calculation showed that the fine structures at the lower energ y region up to E x ~ 16 MeV 


observed in the experimental IS monopole strength function 


33, 


52] surely correspond to 


cluster states and are rather satisfactorily reproduced by its calculation 33], although the 
fine structures are much difficult to be reproduced by the mean-held calculations such as 


RPA and QRPA 54] including the relativistic approach 55] together with beyond mean- 
held calculations |s6] . These results were supported by a recent microscopic calculation with 
an extended 12 C+a cluster model H] based on the antisymmetrized molecular dynamics 
(AMD). On the other hand, the recent experiments on the inelastic form factor to the Hoyle 


state (Oj) as well as its monopole matrix element by the 12 C(e, e') reactions 


HE] 


have 


conhrmed the old experimental data [6i 
by the 3a cluster-model calculations (lOl. 
(FMD) approach 


Their data are known to be well reproduced 


, i 13 ] as well as the fermionic molecular dynamics 

n 

[76| . This fact shows that the a cluster model favorably compares with 
the FMD calculation in describing the monopole transitions, although the model space of 
the 3a cluster model is limited compared with that of FMD. 

The IS monopole transition strengths are also useful to search for cluster states in neutron 
rich nuclei. Quite recently the enhanced monopole strengths in 12 Be, predicted by the gener¬ 


alized two-center cluster model (GTCM) 58 60], have been observed in the breakup-reaction 


experiment using a 12 Be beam 61]. According to their results, the enhanced monopole state 
observed corresponds to the 0 + state at E x = 10.3 MeV with an a+ 8 He cluster structure. 
Thus the IS monopole transition strengths indicate to be a good physical quantity to explore 
cluster states in light nuclei. 

The purpose of the present paper is to study the structure of 13 C from the viewpoint of a 
clustering, in particular, paying our attention to the 1 /2~ states together with the l/2 + ones 
up to around the 3a + n threshold (E x ~ 12 MeV). There exist five 1/2' states (E x = 0.0, 
8.86, 11.08, 12.5 37], and 14.39 MeV) and three l/2 + states (3.09, 11.00, and 12.14 MeV) up 
to E x ~ 12 MeV 62], where several cluster disintegration thresholds are open: 12 C(g.s)+n 
(E x = 4.95 MeV), 12 C(Hoyle)+n (12.60), 9 Be(g.s)+a (10.68), 8 Be+ 5 He (13.2), and 3 a + n 
(13.22) etc. The motivations why we focus on studying the structure of the 1 /2 =b states are 
threefold. The first motivation is originating from the results by the shell model calculations. 
The shell model calculation by Millener et ah 63] and the no-core shell model one by Navrat.il 
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et al. 


67] reasonably reproduced the energy levels of the lowest three 1/2 states and the 


lowest two 1/2 states, respectively. The 1/2 12 states are interpreted as p-shell levels 


dominant configurations are SU(3)[/](A, fi) = [441](0,3) and [432](1,1), respective 


based on the experiments with one-nucleon and two-nucleon pickup reactions 
the Oku shell-model calculation 


64, 


A 


the 


63j, 


65] and 


shell-model state 


66] in 1960’s, while the l/2 3 state is assigned as the 2Tilo 


63]. However, the CO transition matrix elements to the l/2 2 3 states 


measured 


state 


3 


by the 13 C(e,e / ) experiments 68], which are the same order as that of the Hoy 


17|, [62], are much difficult to be reproduced within their shell-model framework 63], 


where CO denotes a longitudinal electric monopole transition. According to their results, the 
difficulty with the transition to the 1/27 8-86 MeV state, in particular, is that the initial and 
final p-shell wave functions have different intrinsic spins, S' = 1/2 and S = 3/2, and thus 
the calculated CO form factor is much smaller than the observed, due to the characteristic 
of the CO transition, i.e. conserving the intrinsic spin. As for the 1/2^ 11.08 MeV state, 
the calculated form factor is also much too small, although its state is described by the 2hoo 
shell-model wave function with dominant intrinsic spin S = 1/2. One of the possibilities to 
solve these difficulties in the shell-model approach may be to take into account the mixture 
of the (0 + 2)Tiuj configurations in each 1/2“ levels including the ground state. However, 
there are no papers reproducing the experimental CO matrix elements with the (0 + 2)huo 
shell-model calculations as far as we know. 

In addition to the experimental CO matrix elements, the IS monopole transition rates 
of 13 C for the lowest three excited 1/2“ states have been reported with the inelastic a 
scattering on the target of 13 C by Kawabata et al. 37]], and their experimental values are 
comparable to the single particle one 8]. In the light of the fact that the monopole strengths 
are one of the good physical quantities to explore cluster states in light nuclei as mentioned 
above, those experimental results suggest that the low-lying excited 1/2“ states have an 
aspect of a-clustering in 13 C. Candidates for the molecular states with K n = 3/2 


and linear-chain configurations of three a clusters in 13 C 69 


i± 


69 


72 


7lJ 


74] have been proposed 


m 


highly excited energy region above the 3a + n threshold. However, the cluster structures 
of the low-lying states including the l/2 ± ones up to around the 3a + n threshold have not 


been studied well, although some preliminary results have been reported 38, [69], [75|. Thus, 
it is important to study the a-clustering aspects in the 1/2” states of 13 C, focusing on the 
CO transitions and IS monopole transitions. 
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The second motivation in the present study is related to the recent progress of the struc¬ 
ture study of 12 C as mentioned above, which has disclosed richness in the structure of 12 C. 
Thus it is very intriguing to study what kinds of structure changes happen in 13 C when 
an extra neutron is added into 12 C, which has the shell-model-like (0/), 3a-gas-like (0^), 


inear-chain-like (O 4 ) structures etc. As men- 


31] are known to describe nicely the structures 


higher-nodal 8 Be(0 + )+a cluster (O 3 ), and 
tioned above, the 3a cluster models [27|, 2!8|, 
of the 0/, O 2 , O 3 ", and O 4 states of 12 C (see also Fig. [2] in this text). On the other hand, the 
third motivation is to explore the Hoylc-analogue states in 13 C, the wave function of which 
is described by product states of constituent clusters such as (0S') 3 (S') n etc. They might 
appear around the 3a + n threshold. 

In the present study, we take a four-body 3a + n OCM with the Gaussian expansion 


method (GEM) 


[77, a 


the model space of which is large enough to cover the 3 a + n 


gas, the 12 C+n cluster, 9 Be+a cluster, and 8 Be+ 5 He cluster configurations as well as the 


shell-model configurations, 
complex scaling method [27 
lying structure up to E x 


etc. Combining OCM and GEM provides a power: 

n n n 

light nuclei [15|, |32j, ]33|J as well as light hypernuclei 


t should be reminded that the 3a OCM calculation with the 
28] succeeded in reproducing well the structure of the low- 
12 MeV including the recently observed O 3 , O 4 , and 2 2 states 

hi method to study the structure of 
7], [79j], because of the following three 
points: 1) the Pauli-blocking effect among the clusters is adequately taken into account by 
OCM, 2) the GEM covers an approximately complete four-body model space 77, Zsj], and 
3) the cluster disintegrated thresholds of the 12 C( 04 " ,2^",4^,0^,3 — , 1 — )-h?z, 9 Be(3/27,l/27)+a, 
8 Be(0 + ,2 + ,4 + )+ 5 He(l/2~,3/2 - ), and 3a+n channels etc. are reasonably reproduced. The 
first point indicates that the OCM-GEM framework can describe the shell-model-like com¬ 
pact structure, for example, in the ground state of 13 C. For the third point, the reproduction 
of the cluster disintegrated thresholds is very important to discuss the cluster structures in 
13 C from the viewpoint of the Ikeda threshold rule 80]. The present framework can explic¬ 
itly treat a strong parity dependence of the a-n potential | 8 l[], demonstrated by the a-n 


scattering phase shifts 82, 83]: the negative-parity potential is attractive enough to make 
resonant states (3/27, 1/27, 7 /2]f and 5/27) of 5 He, wheres the positive-parity potential 
is weakly repulsive. This strong parity dependence plays an important role in producing 
the cluster states of 13 C, as will be discussed later. With this 3a + n OCM framework, we 
investigate the cluster structures and monopole excitations as well as the Hoyle-analogue 
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FIG. 1: Jacobi-coordinate systems in the 3 a + n model. (a)~(d) correspond to the respective 
(p-th) Jacobi coordinate systems (p = 1 ~ 4). See the text. 


states in 13 C. 

The present paper is organized as follows: In Sec. HU the four-body 3a + n OCM is 
formulated. Results and discussions are devoted to Sec. IIII1 Finally we present a summary 
in Sec. |IVl 


II. FORMULATION 

In this section, we present the formulation of the 3a + n four-body OCM. 

A. 3 a + n OCM 

The total wave function of 13 C, Tj( 13 C), with the total angular momentum J and total 
isospin T = | in the 3a+ n OCM framework is expressed by the product of the internal wave 
functions of a clusters 0(a) and the relative wave function T,/( 13 C) among the 3a clusters 
and the extra neutron, 

$j( 13 C) = < 4>j( 13 C)0(a 1 )0(a 2 )0(a 3 ). (1) 
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The relative wave function 4j( 13 C) is expanded in terms of the Gaussian basis as follows: 


$ 


f!3 


c ) = £££/<f> <p) )« , S>(‘' (p) ), 


P= l c (p) 


= S.» 




r (P) 




/» 


M<M id c)> = o, 


( 2 ) 

(3) 

(4) 


where we assign the cluster number, 1, 2, and 3, to the three a clusters (spin 0), and the 
number 4 to the extra neutron (spin |). (z/ p )) denotes the relative wave function with 

respect to the p-th Jacobi-coordinate system of the four-body 3 a + n model (either of coordi¬ 
nate type of K or H) shown in Fig. [1} in which r\ p \ and r-f^ are the Jacobi coordinates. 
S a stands for the symmetrization operator acting on all a particles obeying Bose statistic, 
and (n) is the spin function of the extra neutron. 0(a) denotes the intrinsic wave function 
of the a cluster with the (Os) 4 shell-model configuration. The angular momentum channel 
for the p-th Jacobi-coordinate system is presented as c^ = £^\ t'jf' 1 (p) , \ , 

where £^ p> denotes the relative orbital angular momentum with respect to in Fig. [Tjwith 
L 03 — ^- 2 ' + and stands for the total orbital angular momentum of the 3 a + n 
system with the total intrinsic spin 4. z/G denotes the set of size parameters, u[ p \ and 


of the normalized Gaussian function, <pe(r, u) = N^(u)r e exp(— ur 2 )Yc(r ) 1 and v is taken 
to be of geometrical progression, 


Gi = l/b„, 


bn. - brnlnd 


n —1 


n — 1 ~ n max . 


(5) 


ft is noted that this prescription is found to be very useful in optimizing the ranges with a 
small number of free parameters (6 m i n , a, n max ) with high accuracy 77], 78|. Equation ()4|) 
represents the orthogonality condition that the total wave function ((2D should be orthogonal 
to the Pauli-forbidden states of the 3a + n system, uf' s, which are constructed from the 
Pauli forbidden states between two a particles in OS', 0 D, and IS states together with 


those between a particle and extra neutron n in OS state 


84]. In the present study, the 


Pauli-forbidden states are removed with use of the Kukulin’s method 85] as shown later. 
The ground state with the dominant shcll-model-likc configuration (Os) 4 (Op) 9 , then, can be 
properly described in the present 3a + n OCM framework. 

The 3a + n Hamiltonian for $j( 13 C) is presented as 

4 3 3 

n = ^T,-T cm + £ V 2a (i,j) + ±V an (i, 4) 

i=l i<j=1 i=l 



















( 6 ) 


+ V3a(l, 2, 3) + V2an(*) 4) + V3 an (l, 2, 3, 4) + Vp au li, 

*<J=1 

where Tj, V 2 a (I 4 n ), (Wan), and Wcm stand for the kinetic energy operator for the i- 
th cluster, a-a (a-n) potential, three-body potential among the three a particles (the two 
a particles and extra neutron), and four-body potential, respectively. The center-of-mass 
kinetic energy (T cm ) is subtracted from the Hamiltonian. 

The effective a-a potential V 2a is constructed by the folding procedure from an effective 
two-nucleon force including the proton-proton Coulomb force. Here we take the Schmid- 


Wildermuth (SW) force 


as the effective NN force. It is noted that the folded a-a 


potential reproduces the a-a scattering phase shifts and the energies of the 8 Be ground- 
band state ( J n = 0 + — 2 + — 4 + ) within the framework of the 2 a OCM. As mentioned in 
Sec. HI the 3a OCM_calculations with the complex-scaling method by Kurokawa et al. 0 
and Ohtsubo et al. 28], who use the a-a folding potential based on the SW force together 


with the phenomenological effective three-body force V% a depending on the total angular 
momentum of 12 C, reproduce the 2 £ and 4^ states, which have been recently identified 
above the Hoyle state 2d-25], and the O 3 " and O 4 states, which have quite recently been 


observed in experiments 


22], The origin of V 3 a is thought to be due to the state dependence 


of the effective nucleon-nucleon interaction and an additional Pauli repulsion arising from 
exchanging nucleons among the 3a clusters. In order to incorporate the angular momentum 
dependence of Wa given in Refs. 27,(28] into the present 3 a + n OCM framework as possible, 
we take the following simple one-range Gaussian 


V 3 a(l, 2 > 3 ) = V 3 a exp [-$, a (pf 2 + p % 3 + P 31 )] ( X + V 3 aL 


(7) 


where the operator L 3a stands for the total orbital angular momentum operator of the 3a 
clusters in the 3 a + n system, and p ab denotes the relative coordinate between the a-th a 
particle and 5-th one. In the present study we use [5 = 0.15 fin -2 , = 31.7 MeV, and 
V 3 a = (63.0/31.7 — 1)/6, where the value of /3s a is the same as that used in Refs. [27, 28]. 
These parameter sets give the strength of (l + = 31.7, 63.0, and 136.0 MeV 

for the total angular momentum L 3a = 0, 2, and 4, respectively, which are comparable to 
the strengths of V?, a = 31.7, 63.0, and 150 MeV used in Refs. [27, 2^]. The calculated 
energy spectra of 12 C with the 3a OCM using the folding potential based on the SW force 
and three-body force V 3 a in Eq. (J7|) are shown in Fig. [2l The resonant states above the 3a 
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FIG. 2: (Color online) Energy spectra of 12 C obtained by the 3a OCM calculation with the 
three-body force V 3 a in Eq. 0 , compared with the experimental data (see the text for the details). 


threshold are obtained by applying the complex scaling method (CSM) to the 3a OCM, 
although only the O 3 " resonant state is identified with the more sophisticated method, i.e., 
the method of analytic continuation in the coupling constant combined with the complex 
scaling method (ACCC+CSM) proposed by Ref. [27]. Compared with the experimental 
data, they are reasonably reproduced. It is known that the 3a OCM describes well the 
structure of the 0/] 2 , 2^ 2 , 4/( 2 , 3 _ , and 1“ states in 12 C. 

As for the effective a-n potential V an in Eq. () 6 lh we use the Kanada-Kaneko potential 81 
which reproduces nicely the a-n scattering phase shifts with the odd and even partial waves 
in the low-energy region together with the resonant states (J 71 " = 3/27 and I/ 27 ) of 5 He. This 
local potential, which has the parity-dependent term, is constructed within the microscopic 
a+n cluster model based on the resonating group method (RGM). In the Hamiltonian (151) . 
the three-body force V^n of the short-range Gaussian type is phenomenologically introduced 
so as to fit the ground-state energy of 9 Be measured from the 2a + n threshold within the 
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FIG. 3: (Color online) Low-lying energy levels of 9 Be obtained by the 2 a + n OCM calculation 
together with the experimental data (see the text for the details). 


framework of the 2 a + n OCM. The origin of V 2 an is considered to be the same reason as the 
case of V 3 a mentioned above. In the present paper, we take the following phenomenological 
Gaussian-type potential, 


v 2an\l>iJi 






where p ab denotes the relative coordinate between a-th and 6 -th particles, and (i,j) is either 
(1, 2), (1, 3), or (2, 3). For simplicity we take ft™ = P 3 a — 0.15 fm -2 . The calculated energy 
spectra of 9 Be with the 2 a + n OCM are shown in Fig. [3j The resonant states above the 
2 a + n threshold are obtained by applying the GSM to the 2a + n OCM, although only the 
l/ 2 + state just above the 2 a + n threshold is given with the bound state approximation. 
It is noted that the l/2 + state is considered to be a virtual state 87] and thus the CSM 


88 ]. Compared with the experimental data, 


fails to fix its identification as shown in Ref. 
the low-lying spectra in 9 Be are reasonably reproduced. It is known that the 2a + n OCM 
describes well the cluster structure of 9 Be in the low-lying region. 

As for the four-body force V 3 an in Eq. (J6]), a phenomenological one is introduced so as 
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to reproduce the energies of the ground state ( 1 / 2 ^) with respect to the 3 a + n threshold. 
The origin of V 3 an is similar to that in the case of the three-body forces Vs a and V 2 an as 
mentioned above. It should be short-range, and hence only act in compact configurations. 
In the present paper, we take the following phenomenological four-body potential, 


Q =9 (A,m) L” 


S SU(3) 


(9) 


where < f>^ ; ^g(L 7r ) with the total orbital angular momentum L represents the 
SU(3)[4441](A, /i) wave function with the total harmonic oscillator quanta Q (Q — 9). It is 
noted that the present 3 a + n model space can be classified into the SU(3) bases with the 
irreducible representation, (A,/i)Q, with partition [/] = [4441], and the total wave function 
of 13 C with negative (positive) parity in Eq. (J2]) can be expanded in terms of the <h^ J ^g(L 7r ) 
bases with odd (even) Q values. The (A, / 1 ) = (0, 3) basis with Q = 9 is the unique Pauli 
allowed state for the negative-parity state of 13 C with Q = 9, which is equivalent to the 
shell-model configuration of (0s) 4 (0p) 9 . Thus this SU(3) basis turns out to be the main 
component in the ground state of 13 C. For simplicity, the strength of the four-body poten¬ 
tial, Vo (Q = 9), is fixed so as to reproduce the experimental energy of the ground state of 
13 C with respect to the 3 a + n threshold: V 0 (Q = 9) = 2.0 MeV. The expectation value 
of this four-body potential does not exceed 3 % of that of the corresponding two-body and 
three-body terms, even for the ground state with the most compact structure, i.e. being the 
most sensitive to the potential. This type of potential based on the SU(3) framework is also 
used as the three-body potential in the a + a + t OCM calculation of n B 36]. 

The Pauli-blocking operator Vp au n in Eq. (J 6 ]), which is based on the Kukulin’s method 85], 
is expressed as 


Vpauli ~ hill A \ U f)( U fl ( 10 ) 

A—>-oo 

which rules out the Pauli-forbidden a-a relative states (/ = OS', IS) 0 D) and Pauli-forbidden 
a — n relative state (/ = OS') from the four-body 3a —n wave function. In the present study, 
we take A = 10 4 MeV. 

The equation of motion of 13 C with the 3 a + n OCM is obtained by the variational 
principle, 


(^>j( C) — E \ <f>j( C)) =0, 


( 11 ) 
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where E denotes the eigenenergy of 13 C measured from the 3a + n threshold. The energy E 
and expansion coefficients in the total wave function shown in Eq. (J2J) are determined 
by solving a secular equation derived from Eq. (11111 . 

It is instructive to study single-a-particle orbits and corresponding occupation probabil¬ 
ities in 13 C. We define the single-cluster density matrix for a clusters, respectively, as 

p (a> (r,r') = {<M 13 C)j W|i(rf»-r'»Wr< 0) -r)| |<M 13 C)), (12) 

d 1=1 


where r[ G ^ (i = 1,2,3) represents the coordinate vector of the ith a cluster with respect to 
the center-of-mass coordinate of the 3a + n system. The calculated method of p is given in 


Refs, pj, 


151,189], 


90]. The singlc-a-particle orbits and corresponding occupation probabilities 


are obtained by solving the eigenvalue equation of the single-cluster density matrix, 


J drp^(r,r’)f^(r') = (13) 

where the eigenvalue denotes the occupation probability for the corresponding single¬ 
cluster orbit /^ with the argument of the intrinsic coordinate of an arbitrary a cluster 
measured from the center-of-mass coordinate of 13 C. The spectrum of the occupation prob¬ 
abilities provides important information on the occupancies of the single-a-particle orbit in 
13 C. If the three a particles occupy only an single orbit, the occupation probability for this 
orbit becomes 100 %. On the other hand, the single-particle density matrix for the extra 
neutron in 13 C is also defined as 


p {n \ry) = (<M 13 C)| |5(ri G) -r'))(5(ri G) -r)| |$j( 13 C)), 


(14) 


(G) 

where r\ ; stands for the coordinate vector of the extra neutron with respect to the center- 
of-mass coordinate of the 3a + n system. The single-particle orbits and corresponding 
occupation probabilities of the extra neutron are obtained by diagonalizing the density 
matrix in the same manner as the case of the singlc-a-cluster density matrix (see Eqs. 
and (fT3ji ). 

The root-mean-square (rms) radius of 13 C in the present OCM is given as 
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where R a ~ a (R a - S Be, R n - 12 c) presents the rms distance between a and a (the third a and 
8 Be, extra neutron and 12 C) in 13 C. In Eq. (jT5|) we take into account the finite size effect of 
a clusters, where the experimental rms radius for the a particle is used in y (r 2 ) . 

The overlap amplitudes or reduced width amplitude is useful to see the degree of clustering 
in nucleus. In the present paper, we study the reduced width amplitudes for the 12 C+n and 
9 Be+a channels, respectively, defined as 
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yZ7J( r 94 ) = 


/ 3! 

2hL! 


<5(^4 - r 94 ) 


/ 2 

r 94 


94)^9 ( 9Be )<K«) 


$j( 13 C) , (17) 


where r n (r 94 ) denotes the radial part of the relative coordinate between 12 C and n ( 9 Be and 
a). The wave function </>j c ( 12 C) (0j g ( 9 Be)) of 12 C ( 9 Be) with the total angular momentum 
Jc (J 9 ) is obtained with the 3a OCM (2 a + n OCM). The spectroscopic factor S 2 is defined 


as 


r 00 

S 2 = dr\r x (V(r)] 2 , 
Jo 

where y denotes the overlap amplitude. 


(18) 


B. Monopole transitions, El transitions, and Gamow-Teller transitions 

The CO (longitudinal electric monopole) transition matrix element between the ground 
state (I/27) and n-th excited 1/2“ state by the (e, e') reaction is given as follows: 

M(CO, 1/2- - l/2r) = <i 1/K ( ls C)|f;TtZSi(r i -Jl cm ) ;! |* 1 / 2 -( 13c )>, (19) 

i =1 Z 

where r* (i = 1 ~ 13) are the coordinates of nucleons and R cm = -k 77,'= 1 r i is fi ie c.o.m. co¬ 
ordinate of 13 C. This is the same definition as the CO transition matrix clement, which is 
related with the e + e _ pair creation process, for example, the 6.05 MeV 0^ —y 0^ transition in 
16 0. On the other hand, the isoscalar monopole transition matrix element from the ground 
state (I/27) to the n-th excited state (1/2~) is defined as 

44(IS, 1/2“ - l/2r) = <ii /2 ;;( 13 C)li>, - iM 2 |4 1/2r ( 13 C)}. (20) 

2—1 
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In the present study we calculate the El transition rate between the ground state (l/2j~) and 
first l/2 + state (1/2/) together with the Garnow-Teller transition rate between the ground 
state of 13 C and l/2f (3/2^) state of 13 N, 

B(El) = |(l> 1/2f ( 13 C)||yii^|r i -J{ c , 1 ,|y 1 (r I ^B c »)||i 1/2 -( 13 C))| 2 , (21) 

1=1 Z 

B(GT) = |<* 7/ _ 1/2r , 3/2r ( 13 N)|| yy | I* 1/=r ( 13 C)> I 2 . (22) 

2=1 

The wave functions of the ground state (l/2“) and first 3/2“ states of 13 N, 
^ > j / =i/ 2~,3/2“( 13 N), are obtained by the four-body 3 a + p OCM calculation, in which only 
the Coulomb forces between the extra proton and three a clusters are switched on in the 
present formulation of the 3a + n OCM (see Sec. Ill AD . although the details will be given 
elsewhere. 

III. RESULTS AND DISCUSSION 

Figure H] shows the energy levels of 1/2“ states together with those of l/2 + states in 13 C 
obtained by the four-body 3 a + n OCM calculation. We found that five 1/2“ states and five 
l/2 + states come out as either bound states against particle decays or quasi-bound states. 

A. 1/2“ states 

First we discuss the structures of the Eve 1/2“ states together their isoscalar monopole 
excitations and CO transitions. 

1. Structures of 1/2“ states 

The l/2]~ state, located at E — —12.3 MeV measured from the 3a + n threshold, is the 
ground state of 13 C. Its calculated r.m.s. radius is Rn = 2.4 fin (see Table |T|), the value of 
which is in correspondence with the experimental data (2.46 fm). According to the analysis 
of the wave function, the main component of this state is the SU(3) irreducible representation 
[/]( X, fi) = [4441] (0,3) of the (0s) 4 (0p) 9 configuration with the lowest harmonic oscillator 
quanta Q — 9 and its dominant angular momentum channel is ( L,S)j = (1, |)i, where L 
(S) denotes the total orbital angular momentum (total intrinsic spin): The component of 
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FIG. 4: (Color online) Energy levels of the 1/2 and l/2 + states of 13 C obtained by the 3a + n 
OCM, compared with the experimental data. The experimental levels of the l/2j ~ 2 3 5 and l/2 /~ 2 3 


states are taken from Ref. 


62j ] and from Ref. 3^] for the 1 /2 4 state, respectively. The threshold of 


the 9 Be(5/2 )+a channel at E x = 13.1 MeV, located between the 9 Be(l/2 + )+a and 9 Be(l/2 )+a 
channels, is presented by the dashed arrow in the left hand side of the panel. 


the SU(3) (0, 3) state (0 hu) is as large as 61 %, and the remaining comes from the 2hco state 


(18 %), ATtoj (11 %), (5 %), and higher fiu (5 %). In the cluster mode 


other than the 0 hu one correspond to a-type ground-state correlations 


the components 


33 ], It is noted 


that the SU(3) state, [4441] (0, 3), is the lowest good-spatial-symmetry state in the SU(3) 
model of 13 C. Reflecting the fact that the nuclear force favors a good spatial symmetry, the 
ground state of 13 C has the dominant SU(3)-like nature, although it has significant a-type 
ground-state correlations in the present cluster model. 

The spectroscopic factors and overlap amplitudes of the 12 C+n and 9 Be+a channels 
defined in Eqs. (TT6ll ~ (fl8|l are useful to see the structure of the 1/2R state. Their results are 
shown in Figs. [5](a) and|6](a). The values of the spectroscopic factors and radial behaviors of 
the overlap amplitudes of the 1/2R state can be explained qualitatively by the SU(3) nature 
of the state. The fact that the spectroscopic factors of the 9 Be+a channels are non-zero 
( S 2 ~ 0.2) indicates that the ground state has not only the mean-held degree of freedom 
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TABLE I: Excitation energies (E x ), r.m.s. radii (R), CO transition matrix elements [M(C0)], 
isoscalar monopole transition matrix elements [Ad(IS)] of the excited 1/2“ states in 13 C obtained 
by the _3 a + n OCM calculation. The experimental data are taken from Refs. [62 
Ref. 


68] and from 


37] for the l/2 4 state. The finite size effects of a particle taken into account in estimating R 


with the 3 a + n OCM (see Ref. jl5] for details). 



E x [MeV] 

Experiment 

R [fm] Af(CO) [fm 2 ] 

E x [MeV] 

3 a + n OCM 

R [fm] Af(CO) [fm 2 ] 

Ad (IS) [fm 2 ] 

1/27 

0.00 

2.4628 


0.0 

2.4 



1/27 

8.86 


2.09 ± 0.38 

11.7 

3.0 

4.4 

9.8 

1/27 

11.08 


2.62 ±0.26 

12.1 

3.1 

3.0 

8.3 

1/27 

12.5 


No data 

15.5 

4.0 

1.0 

2.0 

1/27 

14.39 


No data 

16.6 

3.7 

2.0 

3.3 


but also a-cluster degree of freedom. This will be discussed in detail in Sec. IIII A 21 

The ground-state wave function of 13 C can also be studied by calculating the Gamow- 
Teller transition rates 5(GT) between the 13 C ground state (I/ 27 ) and 13 N states 
(l/2]“,3/27) in Eq. (j221) . together with the El transition rate B(E1) between the ground 
state and first l/ 2 + state of 13 C in Eq. (l 2 T]h The calculated 5(GT) values together with 
the experimental data are given as follows: I? cal (GT) = 0.332 vs. T> exp (GT) = 0.207 ± 0.002 
for the transition from the 13 C(l/27) state to 13 N(l/27), 5 cal (GT) = 1.27 vs. I? exp (GT) = 
1.37 ± 0.07 from 13 C(l/27) to 13 N(3/27). The calculated results are in agreement with 
the experimental data within a factor of 1.5. On the other hand, the calculated value of 
B (El : I /27 —> 1/2/) is 2.0 x 10 ~ 3 fm 2 in the present study, while the experimental value is 
14 x 10 -3 fm 2 . This enhanced El transition rate has been pointed out by Millcner et al. 63], 
where the result of the shell model calculation is i?(El) = 9.1 x 10 ~ 3 fm 2 , which is about 
two-third of the experimental value. 

The four excited 1/2“ states, 1/27, 1/27, 1/27, and 1/27, have the larger nuclear radii 
(3.0, 3.1, 4.0 and 3.7 fm, respectively) than that of the ground state (see Table [J). As for 
the I /27 state, the dominant spectroscopic factor is the 9 Be(3/27(g.s))+a channel (see 
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3/2, 5/2, 3/2, 5/2, 3/2, 5/2, 



FIG. 5: (Color online) Spectroscopic factors S 2 of the 12 C(JJ)+n channels (J<5 = 
0+, 1J~, 2+, 3j”, Oj) and 9 Be(Jg)+a (Jg = 3 / 2)j~, 1/2^, 5/2j~, 1 /2^) in the five 1/2“ states of 13 C 
dehned in Eq. (fT 8 l) . 


Fig. E](b)), although one sees non-negligible contributions from the 12 C(3)”)+n channel etc. 
From the analyses of the overlap amplitudes of the 1/2^" state shown in Fig. [HXb) , the relative 
wave function between 9 Be and a in the 9 Be(3/2j“(g.s))+a channel, drawn by the real line, 
has a two-node D -wave radial behavior ( 2D ), where its magnitude is rather small in the 
inside region (r < 2.5 fm), and the magnitude of the maximum peak around r = 4 fm is 
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FIG. 6 : (Color online) Overlap amplitudes of the 12 C+n channels and 9 Be+a channels for the 
five 1/2“ states of 13 C defined in Eqs. (fT6l) and In the panels we present only the overlap 

amplitudes with the S 2 factor larger than 0.2 in Fig. [5j 


about twice larger than that of the amplitude in the inside region (r < 2.5 fm), and the tail is 
extended to the region of r ~ 12 fm. It is noted that its radial behavior is much different from 
the 2 Zl-type harmonic oscillator wave function (h.o.w.f.) with the nucleon size parameter 
b = Jh/Muj (M: nucleon mass), in which the magnitude of the three peaks are almost the 
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same and its h.o.w.f. has no longer tail compared with the real line in Fig. [G^b) . The reason 
why the relative orbital momentum is ^94 = 2 comes from the angular momentum coupling 
between the angular momentum of 9 Be (J g = 3/2 _ ) and 9 Be— a relative orbital angular 
momentum (£ 94 = 2 ), J — Jg + ^W These results indicate that the 1 / 2 ^" state has the main 
configuration of 9 Be(3/2j“)+a, in which the a cluster orbits around the 9 Be(g.s) core with 
2D orbit, that is, a 9 Be+a molecular structure is formed. Since the ground state (l/2(f) 
has a 9 Be(3/2h)+a cluster degree of freedom with one-node D- wave (ID) radial behavior 
(see Fig. | 6 ](a)) , the 1/2^ state can be regarded qualitatively as an excitation of the relative 
motion between 9 Be and a in the ground state, i.e., from ID to 2D. However, one should 
note that the excitation from ID to 2D is not merely the 2 hu excitation in the sense of the 
shell model, because the ground state has significant a-type ground-state correlation and the 
radial behavior of the 9 Be(3/2(f )+a relative wave function in the 1/2^ state qualitatively is 
much different from the 2D-type h.o.w.f., as mentioned above. 

On the other hand, the 1/2 3 state has the dominant spectroscopic factor of the 
9 Be(l/2~)+a channel (5 2 = 0.70), and its overlap amplitude has a maximum peak in 
the outmost (r ~ 4 fm). Thus, this state mainly has the 9 Be(l/2~)+a cluster struc¬ 
ture, where the a particle occupies an 5 orbit with 35-like oscillatory behavior around 
the 9 Be(l/2~) core (see Fig. 13(c)), although the l/2(f state has non-negligible components 
of the 12 C( 02 , 2f, 3(f)+n channels with 5 2 ~ 0.2. The ground state (l/2(f) has also the 
9 Be(l/2f)+a cluster degree of freedom with 25 behavior (see Fig. 13(a)). Thus, the 1 /2^" 
state can be regarded as an excitation of the relative motion between 9 Be(l/2 _ ) and a in 
the ground state, i.e., from 25 to 35. 

As for the 1 / 2 J and 1 / 2 g" states, the analyses of the spectroscopic factors and overlap 
amplitudes indicate that the 1/2J state has a 9 Be(3/2”)+a cluster structure with higher 
nodal behavior (3D), where its overlap amplitude has a maximum peak around r = 7 fm. 
On the other hand, the 1/2J state has a 9 Be(l/2 _ )+a cluster structure with higher nodal 
behavior (45), where its overlap amplitude has a maximum peak around r = 6 fm, although 
one sees non-negligible component of the 12 C( 02 ~)+n channel (5 2 ~ 0 . 2 ). 

Here it is interesting to clarify the reasons why the 9 Be+a cluster states together with 
their higher nodal states come out in the excited 1/2 - states, and the 1/2“ state with the 
dominant configuration of 12 C(Hoyle)+n does not appear in the present study. As mentioned 
in Secs. Hland HTAl the 3a OCM, the model space of which is the subspace of the present 3a+n 
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model space, reproduces well the low-lying structure of 12 C including the O 3 and O 4 states 
etc. , which have been recently observed above the Hoyle state (see Fig. [2]). According to 
the 3a OCM analyses, the Hoyle state is characterized by a dominant 3a-gas-like structure, 
where the three a clusters move gas-likely, only avoiding mutual overlap due to the Pauli 
blocking effect, although the Hoyle state has a non-negligible 8 Be+a correlation. On the 
other hand, the O 3 " and 0^ states of 12 C have the 8 Be(0 + )+a structure with higher nodal 
behavior and linear-chain-like structure with the dominant configuration of 8 Be(2 + )+a, 
respectively. With these facts in mind we will first see what kinds of structures appear 
in the 1/2“ state of 13 C, when an extra neutron is added into the Hoyle state. In this 
case the extra neutron should be P -wave with respect to the center of mass of the 3a 
system, because of the parity conservation and angular momentum coupling. As mentioned 
in Sec. Ill A1 the a-n force has the strong parity dependence: The odd-parity a-n force is 
significantly attractive to produce the P -wave resonant states ( J n = 3/2“ and 1 /2 — ) in 5 He 
(a+n) system, while the even-parity one is weakly repulsive. This attractive P -wave a-n 
force makes a bound state in 9 Be (J n = 3/2~) with respect to the 2a+n threshold, having 
a loosely coupled 2 a+n structure, while the J 71 " = l/ 2 + state appears a quasi-bound state 
just above the three-body threshold (see Fig. [3]). Thus, with the addition of the P -wave 
extra neutron into the Hoyle state, the attractive a-n force reduces the size of the Hoyle 
state with the 3a gas-like structure to reinforce the 8 Be+a correlation in the 3a + n system. 
The extra neutron can be bound with the 8 Be part and then the °Be+a cluster structures 
are produced in excited 13 C states. In fact the 1/2^ and l/2g states have the 9 Be(3/2 _ )+a 
and 9 Be(l/2~)+a cluster structures in the present calculation, where the 9 Be(3/2 _ ) and 
9 Be(l/2 _ ) states are the spin-doublet states each other. 

On the other hand, the 1/2J and l/2g states in the present calculation come out as 
the higher nodal states of 9 Be(3/2 _ )+a and 9 Be(l/2 _ )+a, respectively, and the energy 
difference between the 1/2J and 1/2^ states (1/2^" and 1 /2^") is as small as a few MeV 
(see Fig. [4j). These situations are similar to that in 12 C. In fact the 0-j" state of 12 C has the 
8 Be(0 + )+a structure with higher nodal behavior, which appears by a few MeV above the 
Hoyle state (see Fig. [2]). Thus, the l/2jj 5 states can be regarded as the counterpart of the 
O 3 state of 12 C. 
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2. Monopole excitations 


The calculated CO transition matrix elements, M(C0), defined in Eq. (1I9|) are shown in 
Table [B The M(C0) values of the 1/2^ and l/2 3 states are in agreement with the experi¬ 


mental data within a factor of 2 at most. It is noted that the shell model calculation 


63 ] is 


difficult to reproduce both the CO matrix elements. These facts support that the 1/2-7 and 
1/2^ states have cluster structures. There are no data for the M(C0) values of the l/2^ 5 
states, although their calculated values are smaller than those of the l/ 27 3 states. 

The isoscalar monopole matrix elements, .M(IS), defined in Eq. (l20j) are also shown in 
Tabled] The calculated results of A4(IS) are about twice larger than those of M(C0) for the 
four excited 1 / 22345 states. The reason is given as follows: In the present four-body 3 a+n 
model, the M(C0) in Eq. (TT9j) has a relation with the isoscalar monopole matrix element 
A4(IS) in Eq. (1201) as follows: 


M( CO, 1/2- - 1/2/) 

= i x .M(IS, 1/2; - 1/2/) - lx <4> 1/2 -( 13 C)|(r, 


«cm) 2 |t 1/2r ( I3 C)), (23) 


where r n represents the coordinate of the extra neutron in the 3 a+n model. The contribution 
from the second term in the right hand side in Eq. (f23ll is much smaller than that from the 
first term for the four excited 1/2^ 3j4)5 states in the present calculation. Thus, the relation of 
M( CO) ~ (1/2) xM(IS) is realized. The total value of the energy weighted sum rule (EWSR) 
for the isoscalar monopole transition in 13 C is given as x 13 x R 2 = 6.5 x 10 3 fm 4 MeV, 
where R denotes the nuclear radius of the ground state of 13 C, and m stands for the mass 
of nucleon. Then, the percentages of the energy weight strength to the isoscalar monopole 
EWSR value are, respectively, 13%, 11%, 1%, and 2%, for the l/2^ 34)5 states. The total 
sum of the percentages amounts to be about 25%, the value of which is comparable to the 
case of 16 0 (~ 20%) [8, 33]. 

The preliminary experimental data of the isoscalar monopole matrix elements in 13 C 


have been provided by the inelastic a scattering experiments, performed at RCNP [37], |91]; 
A4 exp (IS) = 6.1 fm 2 ( 6 %), 4.2 fm 2 (3%), and 4.9 fm 2 (5%), for the I/ 27 , l/ 2 3 , and 1 / 2 J 
states, respectively, in which the values in parentheses denote the percentages of the energy 
weight strength to the EWSR value. It is noted that these experimental Al(IS) values are 
comparable to the experimental data of the isoscalar monopole matrix element to the 0 ^ 
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state (Hoyle state) in 12 C, At (IS) = 10.7 fm 2 , together with that to a 8 He+a cluster state 
in 12 Be, At (IS) = 7.1 ± 1.0 fm 2 observed in a recent breakup-reaction experiment using a 
12 Be beam at 29 MeV/nucleon 61]. The sum of the experimental percentages for the 1 / 2^34 
states of 13 C amounts to about 14% of EWSR. The calculated results of the At (IS) values 
together with the sum of the percentages to the isoscalar monopole EWSR value obtained 
by the present 3a + n OCM are in correspondence with the preliminary experimental ones 
within a factor of 2 at most. 

As discussed in Sec. IIII Al the four excited 1 / 2 “ states, 1/2-7, 1/27, 1/2J and 1/27, 
have the following characteristic structures, 9 Be(3/2~)+a, 9 Be(l/2~)+a, 9 Be(3/2~)+a with 
higher nodal behavior, and 9 Be(l/2 _ )+a with higher nodal behavior, respectively, although 
the ground state (1/2R) has a shell-model-like structure. Here it is interesting to investi¬ 
gate the mechanism of why these characteristic 1 / 2^345 states are exited by the monopole 
transitions from the shell-model-like ground state. 

For this purpose, it is instructive to demonstrate that the SU(3) wave function, 
|(Os) 4 (Op) 9 [4441](0, 3)z,=i— <£) %i(n)) J= i-, which is dominant in the ground state of 13 C in 
the present study (see Sec. IIII A II) . is mathematically equivalent to a single-cluster-model 
wave function of 9 Be+a with the total harmonic oscillator quanta Q = 9, 


|(0s) 4 (0p) 9 [4441](0, 3 )l =1 - < 8 > X±( n )) J= i~ 


= AA 


= N 1 


/ 9!4! 
13!" 

I3T 


A \ u 42 (r 94 )0 3 -( 9 Be) 0 (a) 


A \ u 40 (r 9A )(f) j=k - ( 9 Be)0(a) 


J=l 


J=i 


(24) 

(25) 

(26) 


where Ns 1 are the normalization constants, and 0 (a) represents the internal wave function 
of the a cluster with the (Os ) 4 configuration. 0.,( 9 Be) stands for the internal wave function 
of 9 Be with the angular momentum j, 0j( 9 Be) = |(3,1 )e<8>X±( n ))ji the spatial part of which 
belongs to the SU(3) representation (A, ju) = (3,1) of the (0s) 4 (0p ) 5 configuration. The 
relative wave function between the 9 Be and a clusters in Eqs. (|25|) and (|26|) is described by 
the harmonic oscillator wave function Uqlm{ v 94 ) = UQL{r 9 4)Y LM (r 9 4 ) with the node number 
n — (Q — L)/2 and orbital angular momentum L, where r 94 denotes the relative coordinate 
between the 9 Be_and a clusters. One can prove Eqs. (T25li and (1261) with help of the Bayrnan- 


Bohr theorem 


53]. Equations Q25j) and (126]) mean that the ground state of 13 C has 9 Be+a 


cluster degrees of freedom as well as the mean-hied degree of freedom. We call this the dual 
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nature of the ground state [8]. This dual nature is also realized in the ground state of 16 0, 
12 B, and n B etc. 

The operator of the isoscalar monopole transition of 13 C, 0(IS, 13 C) = X)l 3 (rj — R cm ) 2 , 
in Eq. (177j) is decomposed into the internal parts and the relative part, 


£>(IS, 13 C) = C>(IS, 9 Be) + C>(IS, a) + 


9x4 , 
13 r94 


(27) 


where C?(IS, 9 Be) and C?(IS, a) stand for the isoscalar monopole operator of 9 Be and a clus¬ 
ters, respectively. The relative part, (9 x 4/13)r*g 4 , can excite the relative motion between 
the 9 Be and a clusters (2 Tioj excitation). In the case of the isoscalar monopole excitation 
from the ground state described in Eq. (1241) to 9 Be+a cluster states, one can easily prove 
that the monopole matrix elements in Eq. (1201) is originated from only the contribution from 
the relative part, (9 x 4/13)r*g 4 , in Eq. (1771) . This proof is similar to the case of 16 0, in 
which the 0^ 3 states with the 12 C(0^,2 4 ")+Q! cluster structures are excited by the monopole 


transitions from the ground state with doubly closed shell structure 8, l33| . Consequently 
the 9 Be(3/2 - ,l/2 - )+a; cluster degrees of freedom embedded in the ground state are ac¬ 
tivated by the monopole operator and then the l/2^" 3 states with the 9 Be(3/2 _ ,l/2 _ )+a: 
cluster structures are excited by the monopole transitions. In fact one can see that the 
overlap amplitude of the 9 Be(3/2 _ )+a channel in the 1/2^ state shows 2.0-like oscillatory 
behavior, which has one node higher than that in the ground state having ID-like relative 
motion (see Figs. [61(a) and (b)). In addition, in the l/2 3 state, the overlap amplitude of 
the 9 Be(l/2“)+a channel has 3S'-likc oscillatory behavior, while that in the ground state 
has 2S'-like one (see Figs. K a) and (c)). Therefore, the l/2 2 and l/2 3 states, which have 
the dominant 9 Be+a cluster structures, can be excited by the isoscalar monopole operator 
from the SU(3) (0,3) state. However, according to the present calculation, the isoscalar 
monopole matrix elements of the l/2 23 states from the SU(3) (0,3) state in Eq. (124j) . 
(l/2 2 - 3 |0(IS, 13 C)|[(0, 3) i=1 - ®x*(n)]./ =i/ 2 -), are less than one-third or further smaller com¬ 
pared with ones shown in Tabled (l/2^ 3 1C>(IS, 13 C)|l/2(f), in which the wave function of the 
1 /2state is one obtained by the present 3 a + n cluster model. As mentioned in Sec. IIII ATI 
the 1/21" state has significant ex-type ground-state correlation: The component of the 0 hco 
basis is 61 % in the present study, and the remaining comes from the higher Tuu bases. We 
found that the calculated isoscalar monopole matrix elements (and CO ones) in Table Q] dom¬ 
inantly conies from the coherent contribution between the cc-type ground-state correlation 
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in the 1/2 7 state and the l/2^ 3 wave functions having spatially developed 9 Be+a cluster 
structures. This mechanism of the excitation of cluster states from the shell-model-like 
ground state by the monopole transition is common to the cases of 16 0, 12 C, n B, and 12 Be 


etc. which have been discussed in Refs. 


33, 


36 


59 


60]. 


As for the isoscalar monopole transitions to the 1/2J and 1/2)7 states, their matrix 
elements are about half or one third smaller than those to the 1/2)7 and 1/2)7 states, although 
the radii of the former states are lager than those of the latter (see Table [T|) . The reasons 
are given as follows: The structures of the 1/2)7 and 1/2)7 states are characterized by the 
higher nodal behavior in the 9 Be(3/2 _ )-a and 9 Be(l/2 _ )-o: relative motions, respectively, 
i.e. one node higher compared with the 9 Be+a relative motions in the 1/2-7 and 1/2)7 states. 
Reflecting this higher nodal character, the maximum peak in the overlap amplitude of the 
9 Be+ct channels in the 1/2states appears at as large as r = 6 ~ 7 fm, while that in the 
1/2J 3 states (the ground state) is located at r ~ 4 fm (r ~ 2.5 fm), as shown in Fig. [61 The 
relative part of the isoscalar monopole operator in Eq. (l27l) can excite the relative motion 
of the 9 Be+a cluster degree of freedom in the ground state by 2hco with no change of the 
relative orbital angular momentum. In other word this operator can mainly populate the 
9 Be+cc cluster states having one node higher in the relative motion of 9 Be+a than the ground 
state. As mentioned above, the 1/2J 5 states correspond to the 9 Be+a cluster states having 
two nodes higher than the ground state. Thus the 1/2J 5 states are not excited strongly by 
the monopole transition compared with the cases of the 1 / 27)3 states. This result indicates 
that the monopole transition to cluster states with larger radius such as the relevant 1/2J 5 
states is not always stronger. 

The main feature of the 1/27 and 1/2J states having one node higher in the 9 Be-a 
relative motion compared with the case of the 1 / 2-7 an d 1/27 states, respectively, can also 
be verified by the strong isoscalar monopole transitions and strong CO transitions between 
them. The calculated matrix elements are given as follows: A4(IS; 1/27 — 1/27) = 28 fm 2 , 
M (CO; 1 / 22 — 1 / 24 ) = 13 fm 2 ; AT (IS; 1/2 J — 1 / 27 ") = 22 fm 2 , AT (CO; l/2 3 —1/2^") = 8 fm 2 . 
These values are about three times larger than those from the ground state (1/2R) to the 
1 / 2 - 7 3 states (see Table [T|). These results support that the structures of the 1 / 2 J and 1 / 2)7 
states are characterized by the higher nodal behavior in the 9 Be(3/2 - )-a and 9 Be(l/2 - )-a 
relative motions, respectively. 

Here we make a comment on the interpretations of the 1/2-7 8 . 86 -MeV state by the 
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present a cluster model and shell model 63, |66]. According to the shell-model calcula¬ 


tion 


63], the l/2 2 state is interpreted as the p-shcll one with the dominant configuration 
32](1,1), based on the experiments with one-nucleon and two-nucleon 


of SU (3) [/] (A, fi) — [ 


pickup reactions 


64 


and the Ohc o shell-model calculation 


66], although the experimen¬ 


tal CO matrix element of the l/2 2 state as well as the 1/2]]" state is not reproduced at all 
in their calculations, and there are no papers reproducing them with the (0 + 2)Kuj shell 
model calculations as far as we know. On the other hand, the present a cluster model 
for the first time has reasonably reproduced the experimental CO matrix elements (and the 
isoscalar monopole matrix elements obtained by the 13 C(a,cd) reaction). However, the CO 
matrix element of the l/2 2 state is overestimated by a factor of about 2 compared with the 
experimental data, while that of the 1/2J state is only 1.1 time larger (see Table HI). 

In addition to the CO matrix elements as well as the IS monopole matrix elements, another 
interesting experimental information on the l/2^" 3 states is the Ml transition strengths to the 
ground state. According to the shell model calculation 63, 66|, the Ml transition strength 
from the l/2 2 state (8.86-MeV) is about twice larger than the experimental data (0.23 
W.u.), while those from the 11.75-MeV (3/2 2 , T = 1/2) and 15.11-MeV (3/2 — , T = 3/2) are 
reproduced well within a factor of about 1.1. It is reminded that the shell model calculation 
claims that the 8.86-MeV, 11.75-MeV, and 15.11-MeV levels are mainly p-shcll states with 
the [432] symmetry. Only the Ml strength from the l/2 2 state (8.86-MeV) is overestimated 
by a factor of about 2. On the other hand, the experimental Ml strength from the 1/2]" 
state (11.08-MeV) is much weaker (0.036 W.u.), and there are no theoretical calculations 
for it, as far as we know. 

From these experimental data and theoretical analyses of the 1/2“ states, the present 
overestimation of the CO matrix element of the 1/2^ state may indicate that the l/2 2 state 
is the admixture of the 9 Be(3/2 - )+a cluster component and the 0 ku shell-model component 
with the [432] symmetry. It is considered that the CO matrix element of this admixture will 
be reduced compared with one obtained by the present a cluster model, and the a cluster 
component in this admixture is dominantly responsible for the monopole matrix element and 
the Ohio shell-model component with the [432] symmetry is mainly attributed to the one- and 


two-nucleon pickup reactions 63, [66]. This admixture may also reduce the Ml transition 
strength of the l/2 2 state, which is overestimated in the shell model calculation [63, 66] - 


Although further theoretical analyses are needed to solve these problems, it is believed that 
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TABLE II: Excitation energies (E x ) and r.m.s. radii ( R ) of the 1/2+ states in 13 C obtained by 


the 3 a + n OCM calculation. The experimental data are taken from Ref. 


621 ]. The calculated CO 


matrix elements M(C0) and isoscalar monopole matrix elements Ad (IS) from the 1/2+ state to the 
nth excited 1/2+ state are shown for reference. 



Experiment 

E x [MeV] 

E x [MeV] 

3a + n OCM 

R [fm] Af(C0) [fm 2 ] 

M(IS) [fm 2 ] 

1/2+ 

3.089 

3.0 

2.6 



1/2 1 

10.996 

11.7 

3.2 

4.9 

11.0 

1/2 1 

12.14 

12.5 

3.1 

4.1 

8.0 

1/2 1 


14.2 

4.0 

2.9 

4.6 

1/2 1 


14.9 

4.3 

0.6 

2.1 


the present study has indicated the importance of the a-clustering aspects in the l/2 2 state. 


B. l/2 + states 


The calculated energy levels of the l/2 + states are shown in Fig. |4] together with the 
experimental ones. We found the five 1/2+ states in the present calculation, although only 
the three l/2 + states have been observed at the present stage. The calculated r.m.s. radii 
(R), spectroscopic factors (5 2 ), and overlap amplitudes are shown in Table UTl Fig. H and 
Fig. [8l respectively. 

The 1/2+ state, which is located just below the 12 C(0+)+n threshold, has the dominant 
S 2 factor of the 12 C(0]/)+n channel (S' 2 rs-/ 0.65). Due to the fact that this state is bound 
by only 1.9 MeV with respect to the 12 C(0+)+n threshold, the overlap amplitude of the 
12 C(0/)+n channel for the 1/2+ state has a very long tail which is extended up to r ~ 12 fm 
(see Fig. Eta)). Thus, the 1/2+ state has a loosely bound neutron structure, in which the 
extra neutron moves around the 12 C(0+) core with IS 1 orbit. The radius of this state is 
R = 2.6 fm, the value of which is by about 10 % larger than that of the ground state 
(R=2A fm). This slightly larger radius is supported by the experimental analysis with the 


differential cross sections of the inelastic scattering 92]. From Fig. EK a ) one notices that the 
overlap amplitudes of 9 Be(3/2 - )+ct and 9 Be(l/2+)+a in the 1/2]/ state has the 2 P- and 
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3/2, 5/2-j 3/2, 5/2, 3/2, 5/2, 



FIG. 7: (Color online) Spectroscopic factors S 2 of the 12 C(JJ)+n channels (J£ = 

0)*", lj”, 2+, 3J(, Oj") and 9 Be(Jg)+a (Jg = 3/2j~, l/2 4 , 5/2j~, 1/2^) in the five l/2 + states of 13 C 
dehned in Eq. (fTHl) . 


25-like oscillations and their largest peaks are located around r — 3 fm. This fact indicates 
that the 1/2+ state has the 9 Be(3/2 _ )+a and 9 Be(l/2 + )+a cluster degrees of freedom. 

As for the 1/2^3 4 5 states, their radii are larger than the ground state and 1/2+ states 
(see Table ED). According to the analyses of the S 2 factors and overlap amplitudes, the 1/2^ 
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and I /23 states have, respectively, the dominant cluster structures of 9 Be(3/2 _ )+a and 
9 Be(l/2 + )+a, in which the a cluster orbits around the 9 Be(3/2 _ ) and 9 Be(l/2 + ) cores with 
3 P and 3 S states, respectively, although one sees the non-negligible contributions from the 
9 Be(l/2 + , 5/2 _ )+a configurations for the 1/2^ state and 9 Be(3/2~)+a ones for the 1/2^ 
state, and the S 2 factors of the 12 C (Hoyle)+a channels is as small as S 2 = 0.10 ~ 0.15 
for both the l/2 ^ 3 states. The overlap amplitudes of the respective dominant channels 
have the maximum peaks around r = 5 fm. As mentioned above, the 1 / 2 + state has the 
9 Be(3/2“)+a and 9 Be(l/2 + )+a cluster degrees of freedom with the 2 P and 2 S behaviors in 
the 9 Be-a relative motion. Thus, the l/2^~ and l/2g" states are regraded as the monopolc- 
type excitation of the relative motions of the 9 Be(3/2”)+a and 9 Be(l/2 + )+a cluster degrees 
of freedom in the 1/2+ state, respectively, i.e. from 2 P to 3 P and 2 S to 3 S. This situation is 
similar to that in the l/2 ^ 3 states, as discussed in the previous sections. In fact the isoscalar 
monopole matrix elements and CO transition matrix elements from the 1/2+ state to the 
1 / 2^3 states in Table UT1 are comparable to or larger than those from the l/ 2 j~ state to the 
1/233 states with the a cluster structures (see Table U). The l/2 3 state with the dominant 
9 Be(l/2+)+a structure is located just above the 3 a + n threshold. It is noted that the 1/2+ 
state of 9 Be also appears just above the 2 a + n threshold (see Fig. [3]) and this state is pointed 


out to be a virtual state 


871.188 ], Thus, it is interesting to study whether the l/2 3 state of 


13 C turns out to be a virtual state or not with imposing proper boundary conditions for the 
present four-body 3a + n model in near future. 

On the other hand, the 1/2^ state is identified as the 9 Be(3/2 _ )+a cluster state with 
the higher nodal behavior, in which the relative wave function between the 9 Be(3/2~) and 
a clusters has one node more than that in the 1 /2^” state. Reflecting the characteristic of 
its higher nodal state, the radius of this state is R — 4.0 fm, the value of which is larger 
than that of the 1/2;/ state (R = 3.2 fm). 

As for the 1/2^ state, this state is interpreted as a 3a + n gas-like state in 13 C as 
shown below. From Figs. [71(e) and [3(e), one sees that the dominant configuration of the 
l/ 2 )j" state is 12 C(Hoyle)+n ( S 2 rs_/ 0.6) with an S'-wave relative motion, strongly coupled 
with 9 Be(l/2 + )+a ( S 2 ~ 0.5). The characteristic features of this state are that 1) the 
nodal behavior of the overlap amplitudes of the 12 C(Hoyle)+n channel together with the 
9 Be(l/2 + )+a one almost disappears in the inner region (r < 3 fm), indicating that the 
Pauli-blocking effect is significantly reduced, 2) the relative orbital angular momentum in 
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FIG. 8 : (Color online) Overlap amplitudes of the 12 C+n channels and 9 Be+a channels for the 
five l/2 + states of 13 C defined in Eqs. (fT 6 l) and In the panels we present only the overlap 

amplitudes with the S 2 factor larger that 0.2 (see Fig. 0. 


the 12 C(Hoyle) +n channel as well as the 9 Be(l/2 + )+a one is S- wave, and thus all the a —a 
relative motions as well as the a — n ones in this state are dominantly S'-wave, and 3) the 
radius of this state is R = 4.3 fm, the value of which is similar to the calculated result of 
the Hoyle state with the 3a OCM. These results indicate that the 1 / 25 “ state has a 3a + n 
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FIG. 9: Radial behaviors of the dominant 5-wave single-a; orbit (real line) with the occupation 
probability of 48 % and the dominant 5-wave extra-neutron orbit, (S) n , (dashed line) with the 
occupation probability of 65 % in the l/2g” state, which are obtained by diagonalizing the single-a- 
cluster density matrix and single nucleon density matrix of the extra neutron defined in Eqs. (USD 
and (ED, respectively. See the text for the details. 


gas-like structure. 

The nature of the a-particle condensate in the 1/2$ state can be investigated with the 
analysis of the single-a; orbit and occupation probability obtained by diagonalizing the single- 
a cluster density matrix defined in Eq. (1130 . It is reminded that in the Hoyle state, typical 
of a condensate state, the occupation probability of the 05-wave single-a; orbit (zero-node 
5-wave Gaussian type) in its state amounts to be 70 %, and thus the Hoyle state is described 
by the product states of a clusters, (05)^, with the probability of 70 % [3]. We found that 
the dominant 5-wave single-a; wave function in the 1/2$ state is the 05-type (see Fig. ED, and 
the occupation probability of this orbit amounts to be 48 %. The radial behavior shown in 
Fig. [9] is similar to that of the single-a;-cluster 5-wave orbit in the Hoyle state (see Fig. 4(a) 


in Ref. 


15]). The total occupation probability of the 5 wave, i.e., the total sum from 05 


and 15 etc., amounts to be 57 %, whereas that of the P wave stands at 14 %. On the 
other hand, we found that the extra neutron orbit with the largest occupation probability 
obtained by diagonalizing the single-particle density matrix for the extra neutron defined in 
Eq. (1T4D is 5-wave and the occupation probability amounts to be as large as 65 % (hereafter 
we call this orbit (5) n ). The radial behavior of the (5)„ orbit is shown in Fig. [9l where 
its nodal behavior almost disappears in the inner region (r < 3 fm), indicating significantly 
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reduced Pauli-blocking effect. These results indicate that the 1/2^“ state can be regarded 
as an a-particle condensate-like state, (OS ') 3 (S')™, with the probability of 2.09/4 ~ 52 %, 
where the dominator denotes the particle number (three a-clusters and one extra neutron) 
and the numerator is the particle number occupied in the (0S) 3 (S) n configuration, that is, 
3 x 0.48 + 0.65. From Fig. [7] we note that the extra neutron moves at the outer side of 
the 3a-gas-like region, avoiding the Pauli-blocking effect between the extra neutron and 3ct 
clusters. 

The reasons why the occupation probability of the single-a-cluster OS' orbit in the 1/2)/ 
state (48 %) is smaller than that in the Hoyle state (70 %) are given as the follows: The P- 
wave a-n force is attractive to produce a bound state and resonant states in the 2 a+n system. 
This means that the 9 Be(3/2“,l/2~)+a correlations are enhanced in the 3a + n system, as 
discussed in the previous sections (see Sec. IIII Al) . In the 1 / 25 " state the 9 Be(3/2 _ ,l/2 _ )+a 
correlations with P -wave relative motion can be seen, for example, in the calculated results 
of the S 2 factors (see Fig. [71(e)) as well as the occupation probability of the P -wave singlc-a 
orbit being 14 % as mentioned above. This two-body correlations are considered to hinder 
the growth of the a-particle-condensation aspect in the 13 C system. 

Here we present some remarks on positive parity states, in particular, higher spin states. 
The shell model calculation by Millener et al. 63] indicates that the 1/22, 5/2 f, and 3/2^ 
states (including a l/ 2 + state) are dominantly formed by coupling an sd nucleon to 12 C states 
with the [431] spatial symmetry and the former three high-spin states are strongly excited 


by the C3 transitions (see Table II and Fig. 9 in Ref. 63]), although several positive parity 


states are also suggested to be produced by coupling a nucleon to the [44] symmetry states 
of 12 C. The group of the former states can not be addressed in the present 3 a + n cluster 
model. However, an extended cluster model, a + a + ”3IV” + " N" + n , can be addressed to 
study the structure of the high-spin states formed by coupling an sd nucleon to 12 C states 
with the [431] symmetry, where ”3 N” (”1V”) denotes 3 H- 3 He clusters (nucleon), because the 
model space of a + a + ”3IV” + ”N” has the [431] symmetry states of 12 C as well as the [44] 


ones. In the hypernuclear physics, an extended cluster model, a + ”31V” +”1V” 
applied to the study of the structure of 9 Be hypernucleus up to E x ~ 25 MeV 


-A, has been 


a 


93], where 


the high-lying (low-lying) states of 9 Be are dominantly formed by coupling a A particle 
to the [31] ([4]) symmetry states of 8 Be. This hypernuclear extended cluster model has 
succeeded in describing the structures of 9 Be up to E x ~ 25 MeV and also reproducing the 
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excitation functions of the 9 Be(in-flight K , 7 r ), 9 Be(7r + , K + ), and 9 Be(stopped K , 7 r ) 
reactions up to E x ~ 25 MeV, where the single-nucleon parentage to the [31] symmetry 


states (including 
cluster model [7, 


he [4] ones) of 8 Be from the ground state of 9 Be are addressed within the 


93]. Thus, the extended cluster model of 13 C is promising and is one of 


our future subjects. 


IV. SUMMARY 


We have investigated the structure of the l/2 ± states of 13 C up to around E x ~ 16 MeV 
with the full four-body 3 a+n OCM. The model space describes nicely the structure of the 
low-lying states (O^, 2 +, 4^, 3“, and 1 _ ) of 12 C, including the 2 ^, 0j]~, O 4 states above the 
Hoyle state, with the 3ct OCM, together with those of 9 Be, 8 Be, and 5 He with the 2a+n, 
2a, and a+n OCM’s, respectively. We have succeeded in reproducing all the five 1 / 2 “ 
states and three l/2 + states observed up to E x ~ 16 MeV. It was found that the 1/2^ and 
I /27 states have mainly 9 Be(3/2”)+a and 9 Be(l/2 _ )+a cluster structures, respectively, 
while the ground state l/2f has a shell-model-like structure. The 1/2J and 1 /2g" states are 
characterized by the dominant structures of 9 Be(3/2~)+a and 9 Be(l/2 _ )+a with higher 
nodal behaviors, respectively. 

The present calculations for the first time have provided the reasonable agreement with 
the experimental data on the CO matrix elements A4(C0) of the 1/2^ (E x = 8.86 MeV) 
and I /27 (E x = 11.08 MeV) states obtained by the (e, e') reaction, isoscalar monopole 
matrix elements A4(IS) of the 1/2^ (E x = 8.86 MeV), 1/2^ (E x = 11.08 MeV), and 1/2J 
(E x = 12.5 MeV) states by the (a, a') reaction, and r.m.s. radius of the ground state. It is 
noted that the experimental values of A4(C0) and A4(IS) are strong as to be comparable 
to the single particle strengths. The reason why the 9 Be+a cluster states are populated by 
the isoscalar monopole transition and CO transition from the shell-model-like ground state 
has been discussed in detail. We found that this mechanism, which is common to those in 
16 0, 12 C, 11 B, and 12 Be etc., originates from the dual nature of the ground state 0,Q : The 


ground state in light nuclei have in general both the mean-held degree of freedom and 
cluster degree of freedom, the latter of which is activated by the monopole operator and 
then cluster states are excited from the ground state. The present results indicate that 
the a cluster picture is inevitable to understand the low-lying structure of 13 C, and the CO 
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transitions together with the isoscalar monopole transitions are also useful to explore cluster 
states in light nuclei. 

From the analyses of the spectroscopic factors and overlap amplitudes of the 9 Be+a and 
12 C+n channels in the 1/2“ states, dominant 12 C(Hoyle)+n states do not appear in the 
1/2 - states in the present study. This is mainly due to the effect of the enhanced 9 Be+ct 
correlation induced by the attractive odd-wave a-n force: When an extra neutron is added 
into the Hoyle state, the attractive odd-wave a-n force reduces the size of the Hoyle state 
with the 3 a gas-like structure and then 9 Be+a correlation is significantly enhanced in the 
3a+n system. Consequently the 9 Be(3/2~)+a and 9 Be(l/2~)+a states come out as the 
excited states, 1/2^ and l/2g , respectively. On the other hand, higher nodal states of the 
1/22 3 states, in which the 9 Be-a relative wave function has one node higher than that of the 
1/223 states, emerge as the 1/2J and 1/2^ states, respectively, in the present study. It is 
reminded that the 0/ state of 12 C has 8 Be+a structure with higher nodal behavior. Thus, 
the 9 Be+a cluster states with higher nodal behavior, l/2/ 5 , are regarded as the counterpart 
of the 0/ state in 12 C. 

As for the l/2 + states, the 1/2+ state appears as a bound state by 1.9 MeV below the 
12 C(0/)+n threshold. This state dominantly has a loosely bound neutron structure, in which 
the extra neutron moves around the 12 C(0/) core with IS orbit. The calculated radius of 
the 1/2+ state (R = 2.6 fm), slightly lager than that of the ground state (R = 2.4 fm), 
is consistent with the experimental suggestion 92] - It was found that the 1/2/ and 1/2/ 
states have mainly the 9 Be(3/2 - )+a and 9 Be(l/2+)+o: structures, respectively, and their 
radii are around R = 3 fm. These two states are characterized by the strong isoscalar 
monopole excitations from the 1/2/ state. We have discussed in detail the mechanism of 
why the two cluster states are excited by the monopole transitions from the 1/2+ state. On 
the other hand, we found that the 1/2/ and 1/2/ states have dominantly the 9 Be(3/2 - )+a 
structure with higher nodal behavior and 3a+n gas-like structure, respectively, although 
experimentally the two states have not been identified so far. The 1/2/ state with a larger 
radius (R ~ 4 fm) is described by the product states of constituent clusters, having a 
configuration of (0S')/(S') n , with the probability of 52 %. Thus, the 1/2/ state can be 
regarded as an a-condensate-like state. 

It is interesting to study the structure of the higher angular momentum states ( J n = 3/2+, 
5/2+, • • •) of 13 C with the present 3 a -\- n cluster model. In addition it is also important to 
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investigate the structure of the 13 N mirror nucleus with 3 a + p cluster model. The Coulomb 
energy shifts between 13 C and 13 N as well as the decay widths of their excited states and 
Ml transition strengths etc. will offer another interesting information on their structures. 
These theoretical studies are now in progress. The results will be given elsewhere. 
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